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We develop a temperature measurement of an atomic cloud based on the temporal correlations of 
fluorescence photons evanescently coupled into an optical nanofiber. We measure the temporal width 
of the intensity-intensity correlation function due to atomic transit time and use it to determine the 
most probable atomic velocity, hence the temperature. This technique agrees well with standard 
time-of-flight temperature measurements. We conhrm our results with trajectory simulations. 


I. INTRODUCTION 

Optical nanofibers (ONF) have emerged as a nonin- 
vasive probe for spectroscopy, sensing, and cold atom 
physics In the case of cold atomic gases, the sub- 

/im diameter of the ONF allows for insertion into the 
atomic cloud with minimal disturbance to the trapping 
beams. Moreover, the small mode area of the evanes¬ 
cent field around the ONF waist leads to strong coupling 
between the ONF guided mode and atoms near its sur¬ 
face [H]. 

Optical nanofibers have been used to study atomic 
spectra near surfaces [ilEllH] and magneto-optical trap 
(MOT) size, lifetime and temperature (by completely 
different methods than reported here) [101 E] • Having 
a good temperature diagnostic is important, for exam¬ 
ple, in optimizing nanofiber trapping from a cold ther¬ 
mal cloud [IMl]. Standard techniques for measuring 
MOT temperature, such as time-of-flight (TOF) absorp¬ 
tion imaging or dithering the magnetic field gradient m, 
are effective, but constrained experimental environments 
can prevent their use. 

Hybrid quantum systems composed of superconduc¬ 
tors and neutral atoms require cryogenic environments 
that are incompatible with many of the diagnostic tools 
used in laser cooling and trapping. To successfully inter¬ 
face cold atoms and superconducting qubits [I7H23] it is 
necessary to develop tests that do not perturb the cryo¬ 
genic environment. Superconductors are perilously sen¬ 
sitive to optical power and DC and AC magnetic fields. 
The standard temperature measurements of the atomic 
cloud mentioned above are inaccessible in these setups 
due to a lack of optical access or sensitivity to changing 
magnetic fields. Here we present, as part of our program 
to magnetically couple atoms trapped around an optical 
nanofiber to a superconducting resonator [TS] , a new way 
to monitor the temperature of a cloud of cold atoms near 
a nanofiber, using the correlations of fluorescent photons 
emitted into the nanofiber guided mode. 

When the emitters are not stationary, the intensity- 
intensity correlation function depends on their motion 


* rolston@umd.edu 


as well as the geometry of the mode into which they 
emit im [^. The intensity-intensity correlation func¬ 
tion, measures correlations in the fluctuations of 

light intensity, e.g. the photon statistics [55], and can 
reveal both classical and quantum aspects of the light 
and its sources. Among the classical effects is the transit 
time of atoms through an optical mode. Systems such as 
atomic beams [57] [55] , single atoms in a MOT [^ |3D| , 
and a single trapped ion [3T| were used to measure these 
transit-time effects. While bunched and antibunched 
photon statistics have been observed in the light emit¬ 
ted into the ONF guided mode [55H55] . the correlations 
related to atomic trajectories near the ONF have been 
studied only tangentially [33|. Intensity correlations de¬ 
cay with a characteristic time that depends on atomic 
transit. We measure this time for different atomic tem¬ 
peratures. Its dependence on temperature allows for a 
simple model to infer the MOT temperature directly from 
the correlations. 

This paper is organized as follows. Section |TT| out¬ 
lines the nanofiber mode structure, potential shifts, and 
the coupling efficiency of fluorescence photons into the 
nanofiber. Section in provides a general overview of 
intensity-intensity correlations and briefly discusses the 
theoretical considerations for calculating and simulating 
them. Finally Section |IV| presents the experimental re¬ 
sults and compares them to simulations. 


II. THE SYSTEM 


The experiment relies on two main parts: a source 
of cold ®^Rb atoms and an ONF. A MOT provides a 
constant source of slowly moving atoms whose fluorescent 
light can couple evanescently into the guided mode of the 
optical nanofiber. The nanofiber collects the light from 
the atoms and also modifies the local potential for the 
atoms, which move with typical velocities on the order of 
10 cm • s“^. 
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A. Nanofiber mode structure 

Our single-mode nanofiber is a fiber pulled to a small 
enough diameter that all higher-order modes are cut off. 
The mode (HEn) of such an ONF has an intensity profile 
outside of the fiber given by m 

\'E{r)\'^= [K^{qr)+ uKf{qr)+wK^{qr)] , (1) 

where is proportional to the intensity at the fiber sur¬ 
face; Ki is the modified Bessel function of the second 
kind of order i; u and w are constants obtained from 
Maxwell’s equations; r is the distance from the center of 
the fiber; and q = is the transverse component 

of the wavevector, where /3 is the field propagation con¬ 
stant in the nanofiber, and k = 27 :jX is the free-space 
wavevector. The parameter q describes the decay of the 
field in the radial direction. 


B. Atom-surface potential 


We approximate the nanofiber as an infinite dielectric 
plane when calculating the van der Waals potential [SBl- 
155] . so that C/vdw(?’) = C 3 {r — rg)”^ for r > rg, where 
rg is the fiber radius. The Cg coefficient is equal to 
4.94 X lO--*® J • m-3 and 7.05 x IQ-^® J • m-^ for the 
55 'i/ 2 and 5 P 3/2 levels, respectively, of ®^Rb near fused 
silica [HEnHlI]- The infinite-plane approximation is ac¬ 
curate to within 20% for atom-fiber distances less than 
200 nm [12], a distance comparable to the dec ay len gth 
of the evanescent field {q~^ « 188 nm, see Sec. 1VB|. 


To include the effect of retardation, which causes the 
atom-surface interaction to scale as (r —ro)“^ for r 3> rg, 
we use a phenomenological model for the potential that 
smoothly connects the non-retarded (van der Waals) and 
retarded (Casimir-Polder) regimes [71H2]: 


U{r) = 


C 4 

(r - ro)3 ((r - rg) -k Ci/Cz) 


( 2 ) 


where is the Casimir-Polder coefficient. This coef¬ 
ficient is equal to 4.47 x 10“®® J • m“^ and 12.2 x 10“®® J • 
m“^ for the 5S'i/2 and 5 P 3/2 levels, respectively, of ®'^Rb 
near fused silica gniiis] 


C. Potential shifts 


The atom-surface potential shifts the atomic levels de¬ 
pendent on position. The shifts produce a spatially- 
varying absorption (emission) rate [44]: 


Pahs (r) 


P s 

^l + s + 4(^^^±^)" ’ 


(3) 


where r is the position of the atom, s = I /Jsat is the 
saturation parameter (/gat = 3.58 mW • cm“® for a uni¬ 
form sublevel population distribution [IS]), 5 = wl — wg 


is the detuning of the driving (i.e. MOT) beams from 
atomic resonance, and doj (r) = {Ueir) — Ug{r)) jh is the 
atom-surface shift assuming a two-level atom. 

Note that we neglect effects due to interference of the 
MOT beams with each other or due to scattering of the 
MOT beams off the nanofiber in the near field. This is 
justified because the fiber is long relative to these effects 
in one direction so that the atoms see a landscape that 
is, on average, uniform. 


D. Coupling efficiency 

The coupling efficiency of an atom to the ONF is the 
rate of spontaneous emission that couples into the one¬ 
dimensional mode of the fiber divided by the total spon¬ 
taneous emission rate [S] |35] , 

/3(r) = PiD (r)/Ttot (r) . (4) 

Fermi’s golden rule determines the form of /3 (r), which 
follows the spatial variation of Fq. [^ 

Photon detection in the experiment is a joint process 
of absorbing a photon from the MOT beams and emit¬ 
ting it into the nanofiber mode, which is mathematically 
described by the product of the photon emission rate 
in Eq. and the coupling efficiency in Eq. It is the 
position-dependence of this joint probability that allows 
us to obtain information about the atomic motion. 


III. CORRELATIONS 

The intensity-intensity correlation function 

„(.)(,), w) at+ 0) ,5) 

measures the conditional probability of detecting a pho¬ 
ton at a time t from recording another photon. Here 
(•) denotes an average over all t, and, in this discus¬ 
sion, /(t) is the photocurrent or, equivalently, the photon 
counting rate at time t. At its core, (t) characterizes 
the fluctuations in I{t). When measuring fluorescence 
from an atomic ensemble, the function contains contri¬ 
butions from different sources of fluctuations including 
single-atom field-field correlations, single-atom intensity- 
intensity correlations, different-atom field-field correla¬ 
tions, and different-atom intensity-intensity correlations. 
Neglecting correlations between different atoms and as¬ 
suming that they are motionless, we can write (r) 
as [23] 



where N is the average atom number in a particular 
time window, and (/^^(r) and g^^\T) are the single¬ 
atom intensity-intensity and field-field correlations, re¬ 
spectively. For small atom number N, we can observe 

/o\ 

the “antibunching term” g\ {t). 










3 


Laser-cooled atoms are not stationary emitters. The 
resonance fluorescence emitted into the fundamental 
mode exhibits correlations due to transit-time effects re¬ 
lated to the geometry of that mode. The atoms act as 
beacons signaling their position while passing near the 
nanofiber. Accounting for the motion of atoms amounts 
to adding a temporal envelope /(r) to Eq. [27], 

g(2)(r) = 1 -h f{T)\g^j\T)\^ + . (7) 

The function /(r) generally depends on the environment 
and how the emitted light couples to the detection appa¬ 
ratus - it is the shape of this temporal envelope that will 
allow us to extract information about the trajectories of 
atoms moving near an ONE. 

We can relate the width of the correlation function 
(t) to the temperature of the atomic cloud by noting 
that the temperature determines the velocity distribu¬ 
tion of the atoms and the speed of the atoms determines 
the timescale of the interaction with the nanofiber. The 
ONE mode described by Eq. possesses a characteristic 
length scale of 1/q. Dividing this length by a characteris¬ 
tic speed of a Maxwell-Boltzmann distribution of atoms 
at a temperature T, which we take to be the most prob¬ 
able speed, Vp = ^j2kBT/m^ yields a simple relationship 
between transit time and temperature: 


a m 
^ qV 2kBT 


( 8 ) 


where a is an overall scale factor based on the geometry 
of the problem and on our choice of characteristic speed. 
We are not able to find an analytical form for a from 
simple physical considerations, but used simulations to 
understand how geometric details affect a (see Sec. IVD |. 


IV. EXPERIMENT AND RESULTS 
A. Apparatus 

We load the MOT from the low-velocity tail of a back¬ 
ground vapor of ®^Rb atoms produced by a dispenser (see 
details in Ref. m)- We change the intensity and detun¬ 
ing of the cooling beams in order to controllably vary the 
temperature of the atomic cloud between ^ 160—840 fiK, 
as measured by time-of-flight expansion through fluores¬ 
cence imaging. Our ability to determine the atomic cloud 
temperatures is limited by the time-of-flight (TOE) imag¬ 
ing system in our setup combined with the low atom 
numbers for colder MOTs in steady state. While we col¬ 
lected correlation data for ostensibly colder MOTs, we 
only present data for temperatures for which we could 
provide calibration to a known technique. This is an in¬ 
dication that in certain circumstances the signal-to-noise 
ratio of the correlation measurement technique can be 
better than that of TOE. 

The optical nanofiber (ONE) is produced via the flame 
brushing technique [121 115] ■ A hydrogen-oxygen flame 


acts as a local heat source to soften 125 ^m-dlameter, 
single-mode fiber (Fibercore SM800) whose ends are 
pulled with computer-controlled linear motors. This 
method reliably produces fibers of subwavelength diam¬ 
eters with transmission of the fundamental mode above 
99% and as high as 99.95%, allowing them to sustain 
powers of hundreds of milliwatts in high vacuum [48] . 
Based on our fiber-pulling reproducibility, we know the 
transmission is greater than 95%. Relying on repeated, 
destructive measurements of the nanofiber diameter us¬ 
ing a scanning electron microscope (SEM), we estimate 
the diameter of our ONE to be 530 ± 50 nm, with a 1% 
uniformity over a length of 7 mm. This fiber diame¬ 
ter with the stated uncertainty accepts only one guided 
mode, described by Eq. above, at the experimentally 
relevant wavelength of 780 nm. 

This same fiber has been in our apparatus for over two 
years with no noticeable degradation in quality. Rubid¬ 
ium atoms can coat the fiber surface and reduce trans¬ 
mission under operating pressures, but application of a 
thru-fiber heating beam with a power of more than a few 
/rW is sufficient to desorb the atoms within a few ms. 

We glue (EPO-TEK OG116-31) the fiber to a tita¬ 
nium U-shaped mount for stability, and attach the mount 
to a UHV-compatible manipulator system (VG Scienta 
Transax). The manipulator consists of a motorized step¬ 
per motor along one axis and 2D manual translations 
stages along the other axes. This manipulator works in 
conjunction with three pairs of magnetic shim coils to 
optimally overlap the nanofiber waist with the region of 
highest atomic density in the cloud. 

Light that couples into the guided mode is filtered 
at the output of the fiber by a volume Bragg grating 
(VBG, OptiGrate BP-785), a narrow-line interference fil¬ 
ter (Semrock LLOl-780-12.5), and a long-pass color filter 
(Thorlabs FGL645) before being sent to the two fiber- 
coupled single-photon counting modules (SPCMs) (see 
Fig. 0. A field-programmable gate array (FPGA) 
stores and time-tags photon output TTL pulses from the 
SPGMs, which are then post-processed and correlated. 
An Internal clock of 48 MHz sets the minimal time res¬ 
olution to 20.83 ns. The use of two SPCMs circumvents 
problems near zero time delay related to detector dead 
time, typically 50 ns. 


B. Data and fitting 

For this experiment, the MOT beams are on continu¬ 
ously during data acquisition and drive fluorescent tran¬ 
sitions in the atoms. We collect ^ 2.5 x 10^ photon counts 
for each experimental run, corresponding to about 45 min 
of averaging per data point. Time-of-flight imaging mea¬ 
sures the temperature of the atomic cloud before and 
after photon collection. 

Our data is a list of times corresponding to photon de¬ 
tection events, which we use to find g^"^^ (r). We do not do 
any further binning of the data, so that the timing resolu- 






4 



FIG. 1. (Color online) Experimental schematic. A MOT is 
spatially-overlapped with a nanofiber, and the MOT beams 
drive resonance fluorescence that couples into the guided 
mode. This signal is hltered by a volume Bragg grating 
(VBG), bandpass (BP) hlter, and long-pass filter before being 
split by a 50:50 beamsplitter (BS) and sent to two SPCMs. 
TTL pulses from the SPCMs are time-tagged by an FPGA 
and correlated in software. 


tion of 20.83 ns is set by the internal clock in the FPGA. 
While this time resolution obscures details on atomic 
spontaneous emission timescales (tens of nanoseconds), 
it provides good resolution on the timescale of a few mi¬ 
croseconds where the atomic trajectories produce signa¬ 
tures in the correlation function. Measurements using an 
oscilloscope (Tektronix DPO 7054) with finer time reso¬ 
lution allowed us to observe antibunching for low atom 
number. 

Varying the Rb dispenser current allows us to change 
the number of atoms in the MOT, so that we can 
change the average number of atoms interacting with the 
nanohber mode. Fig. shows the transition from an¬ 
tibunched (positive slope after r = 0, estimated atom 
number is ^ 1.4) to bunched (negative slope after r = 0, 
estimated atom number is ^ 6) correlations as we in¬ 
crease the number of atoms fluorescing into the mode of 
the ONF. The timescale of the bunched or antibunched 
feature is set by the internal degrees of freedom of the 
atom and is much shorter than the temporal envelope 
due to atomic motion, which we discuss next. Similar 
results were observed in Ref. [34| . 

Figure [^displays an example of (r) extracted from 
data for an atom temperature of 460 p,K. Note the very 
different timescale than in Fig. so that the atom- 
number-dependent peak or dip corresponds to only one 
data point in Fig. Because this point at zero time 
delay is the only one that depends on atom number, we 
neglect it when fitting the data. The signal has a char¬ 
acteristic width due to transit-time effects, which is the 
result of a position-dependent atom-fiber coupling effi¬ 
ciency combined with moving atoms. An atom at a par¬ 
ticular location will emit into the mode with probability 
proportional to its intensity at that position, and averag¬ 
ing over many atomic trajectories will sample the entire 
mode. In this way, the autocorrelation function contains 
information about the mode in question (the shape of 
(r)) and about the motion of the atoms (the decay 
time of (t))- 

We make a series of approximations to the model of the 



FIG. 2. (Golor online) Second-order correlation function 
g^^\T) for light scattered into the fiber as a function of delay 
time r. The curves show data for low (blue) and high (red) 
Rb dispenser currents, illustrating antibunching and bunch¬ 
ing, respectively. 
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FIG. 3. (Golor online) Second-order correlation function 
(r) as a function of delay time r for an atom temper¬ 
ature of 460 pK. The data (blue dots) are fit (solid red line) 
to 1 -I- /(r) using Eq. |10[ with the residuals displayed in the 
lower plot. 

mode structure before comparing the observed transit¬ 
time broadening to theory. The factors u and w in Eq. 
are small for a fiber radius of 265 nm and wavelength of 
780.24 nm (calculated to be 0.166 and 0.00875, respec¬ 
tively), so we neglect them and keep only the first term, 
which is proportional to Kq. As a further simplifying ap¬ 
proximation we also take the asymptotic form of Kq [51j , 



which is valid in our case. This yields an intensity around 
the nanohber proportional to eyip[—2qr]/2gr. Dehning an 
effective index of refraction, Ueff = /3/A:, we can rewrite 
the propagation constan t so tha t the radial decay pa¬ 
rameter becomes q = Ay/n^j^ — 1, which is 0.66A for our 
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nanofiber. We recast the spatial dependence of the in¬ 
tensity into the temporal envelope in Eq. [271 EH] : 


g-2(|T|/To-|-0.66 k rg) 

/('^) ^ (|r|/ro-I-0.66fcro) ’ 


( 10 ) 


where ro = 265 nm is the fiber radius, A is a fitting pa¬ 
rameter for the overall amplitude, and the absolute value 
reflects the time-symmetric nature of the autocorrelation 
function for stationary processes. The parameter tq rep¬ 
resents a characteristic correlation time (see Eq. 1^. 

The red curve in Fig. ^ shows the best fit to g^ir) = 
1 -|- /(r) because g\ {t) is flat in our experiment over 
timescales longer than the atomic lifetime. For a mea¬ 
sured atomic temperature of 460 /iK, the fit achieves a 
reduced of 1.02 (and a range of approximately 1 — 1.5 
across all datasets). We note that using Eq. 10 for the 
temporal envelope /(t) results in statistically better fits 
than an exponential or Gaussian decay. 

The overall height of the bunched peak depends on 
the absolute knowledge of the background at very long 
time (seconds), which we know can depend on mechani¬ 
cal vibrations of the fiber. The environment acoustically 
and thermally drives these vibrations so that we do not 
have an exact measure of unity. This, combined with the 
signal-to-background ratio of the photon counting rates, 
can explain why the amplitude A does not reach the ex¬ 
pected value of 2 for chaotic light from independent emit¬ 
ters. 


C. To VS. temperature 

We extract best-fit values for tq at different MOT 
atomic temperatures, with each temperature also mea¬ 
sured by standard TOE imaging. Fig. [^ shows a plot 
of the resulting best-fit values tq, where the vertical er¬ 
ror bars are the standard errors from the fit, and the 
horizontal error bars originate from uncertainty in the 
knowledge of the magnification of the imaging system. 
The average atom number for all points in Fig. [^ falls in 
the same range as for the data presented in Fig. i.e. 
1 < IV < 10. The purple line in Fig. is a fit to Eq. 
with the single fit parameter a, and the shaded area rep¬ 
resents the 95% (2cr) confidence bands considering both 
the vertical and horizontal error bars. We observe good 
agreement between the model and the data; the fit has 
a reduced of 1-67, and the overall scale parameter is 
a = 1.46 ± 0.04. The exact value of this scale factor is 
discussed further in Sec. IlVDl 


D. Simulations 


To better understand the physical situation, we per¬ 
form simulations of atomic trajectories subject to New¬ 
ton’s equations of motion |S] . These simulations include 



FIG. 4. (Color online) Extracted tq vs. temperature T, mea¬ 
sured via TOF. The blue circles are experimental data. The 
vertical error bars indicate standard error in the ht of Eq. |10| 
and the horizontal error bars arise from systematic uncer¬ 
tainty in the magnification of the imaging system. The purple 
line is a ht to Eq. [^ and the shaded region is the 95% (2 ct) 
conhdence bands considering both the vertical and horizontal 
error bars. The reduced is 1.67. The red open squares are 
the results of the trajectory simulation, with a single scale 
parameter of 0.77. 


the atom-surface potential and its resultant shift dis¬ 
cussed in Sections IIIBI and III Cl The classical nature of 
the simulations is justified because the smallest angu¬ 
lar momenta present in the system are still ~ 100 times 
larger than h. 

The atoms are started at a radial distance r = 1500 
nm away from the fiber surface. At this distance, the 
coupling is weak due to the rapid decay of the mode 
with length scale 1/q. The axial and radial symmetry 
of the problem allows us to restrict trajectories to the 
x-y plane with initial velocities pointing in one quadrant. 
We sample the speeds from a 3D Maxwell-Boltzmann dis¬ 
tribution before projecting onto this plane. Trajectories 
evolve for either 50 gs or when the atom strikes the fiber 
surface, whichever happens first. 

The coupling efficiency in Eq.j^is a fit to the complete 
solution for a two-level atom [IS]. We also assume that 
the orientation of the atomic dipoles relative to the fiber 
surface is random, so that the coupling efficiency is an 
effective ensemble average. Independent measurements 
confirm that this assumption of random orientations is 
valid for our MOT. 

Photon scattering events are infrequent on microsec¬ 
ond timescales, and when they do occur their effect on 
atomic velocity is negligible. As a result, we assume bal¬ 
listic trajectories. The position-dependent coupling ef¬ 
ficiency in Eq. [^ and the position-dependent emission 
rate in Eq. are calculated at each instant of time 
along a trajectory and multiplied together. This yields 
a time-dependent detection probability for each trajec¬ 
tory. Time-correlating the detection probability of a tra¬ 
jectory with itself produces a signal proportional to the 
intensity-intensity correlation for a single atom. We dis- 
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FIG. 5. (Color online) Simulated correlation time tq vs. sam¬ 
pling angle range A9 for an atomic temperature of 90 fiK. 
The dashed line indicates the critical angle arctan(53/300) in 
the simulation at which not all atoms hit the fiber. 


cretize these time-dependent probabilities onto a mesh of 
50-ns resolution so that calculating the correlation func¬ 
tion becomes a simple array operation. Experimentally- 
measured values for atom temperature are fed into the 
simulation, which is averaged over 5,000 randomly sam¬ 
pled speeds and directions. The resulting correlation 
function is fit to Eq. [I^ in order to extract the decay 
time Tq. 

We first utilize the simulations to address the scale fac¬ 
tor in Eq.[^ Fig. [^displays the dependence of the transit 
time on the angular spread of the atomic trajectories for 
a distribution with temperature 90 fiK. For an atomic 
beam aimed directly at the fiber, we extract a transit 
time of 1.49 /rs, which matches well the calculated time of 
1.43 ns using Eq.j^with a = 1. The transit time increases 
slowly as we increase the angular distribution of trajecto¬ 
ries, until it hits a critical value of arctan(53/300), corre¬ 
sponding to the point after which not all paths intersect 
with the nanofiber. Beyond this angle, atoms then inter¬ 
act with the fiber over distances longer than l/g, and the 
transit time consequently increases further. Fig. [^illus¬ 
trates that the simulation fully samples the interaction 
region with an angular spread of at least tt/G to get rea¬ 
sonable results. Moreover, we note that the ratio of the 
transit time for the fully-sampled simulation to the effec¬ 
tive one-dimensional simulation with no angular spread 


is 1.7. These results suggest that our observed scale pa¬ 
rameter of o = 1.46 ± 0.04 is partially due to angular 
spread in the trajectories. 

We performed simulations for the same temperatures 
measured in the experiment and one additional atomic 
temperature of 100 ^K. The red open squares in Fig. [^ 
display transit times extracted from simulated data that 
were fit to the temporal envelope function, Eq. The 
simulated data are multiplied by a single, overall scale 
parameter equal to 0.77 in order to minimize the least- 
squares distance to the experimental fit. The simulations 
follow the expected trend with temperature and differ 
only by a scale factor of order unity. The discrepancies 
might be explained by the various simplifications made in 
our model, which neglects, for instance, the stochastic na¬ 
ture of the photon absorption/emission process. We note, 
however, that the difference between the experimental 
and simulated data is comparable to other temperature 
measurement methods using optical nanofibers m- 

V. CONCLUSIONS 

We have presented a technique to measure the tem¬ 
perature of a laser-cooled atomic cloud that is applica¬ 
ble to experiments with restrictive environments, such as 
hybrid quantum systems using superconducting circuits. 
The method uses the intensity-intensity correlation func¬ 
tion to extract motion of atoms as they pass through the 
ONF mode and is easily extendable to other photonic de¬ 
vices with different optical mode geometries. This tech¬ 
nique allows mapping of mode structures, which could be 
useful when using the next family of higher-order modes 
to trap atoms near an optical nanofiber |52H56) . 
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